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Summary  of  Research 


The  primary  purpose  of  this  research  was  to  develop  the 
phenomenon  of  extended  x-ray  absorption  fine  structure  (EXAFS)  as  a 
technique  for  investigating  the  structure  of  non-crystalline  materials, 
particularly  amorphous  semiconductors.  A  theory  had  been  developed 
(Sayers,  Lytle,  Stem,  "Advances  in  X-ray  Analysis",  Vol.  13,  p.  248 
(Plenum,  New  York,  1970)  which  treated  EXAFS  as  a  scattering  of  the 
ejected  photoelectron  by  the  near  neighbors  of  the  absorbing  atom. 

Soon  after  the  research  began  a  major  breakthrough  was  achieved  in 
which  the  theoretical  EXAFS  expression  and  experimental  data  could  be 
formally  inverted  into  a  radial  structure  function  using  the  Fourier 
integral  theorem  The  technique  was  further  developed  and  applied 
to  an  investigation  of  the  structures  of  Ge,  Ge02,  and  GeSe^  to  make 
the  important  point  that  the  environment  of  all  the  different  atoms  in 
a  multicomponent  material  could  be  determined.  This  work  was  recognized 
as  a  significant  contribution  to  x-ray  physics  when  Dale  Sayers  was 
awarded  the  1972  Sidhu  award  by  the  Pittsburgh  Diffraction  Conference. 

Dr.  Sayers'  Th.D.  Lhesis,  w*„  ch  described  this  work  in  detail  was  sub¬ 
mitted  as  a  special  technical  report. 

The  surface-active  character  of  evaporated  amorphous  germanium  (E 
was  emphasized  in  a  calorimetric,  diffraction  and  microscopic  evaluation^) 
which  showed  that  3  wt.%  1^0  could  be  adsorbed  when  E  Ge  was  exposed  to 
air  and  could  not  all  (50%)  be  removed  unless  the  sample  was  heated 
above  the  re crys tallization  temperature.  During  recrystallization  electron 
diffraction  patterns  showed  transient  (222)  reflections  (normally  forbidden) 
which  were  interpreted  as  evidence  of  wurtzi  te-type  bonding  in  the 
original  material. 

The  EXAFS  work  on  amorphous  semiconductors  was  continued  on 
GeSe,  GeSe2,  As2Se^,  As2Sj,  and  As2Te^**)  and  reported  at  the  1973 
Conference  on  Amorphous  Semiconductors  at  Garmis  ch-Fartenkirchen . 

This  work  developed  a  normalization  technique  which  allowed  determination 
of  absolute  coordination  numbers  as  well  as  identifying  the  kinds  of 
atoms  in  the  first  coordination  sphere.  For  example,  in  ancrphous 


As^Se^  we  found  the  As  atom  to  have  3  nearest  neighbors  (2Se  and  1  As) 
while  Se  had  2  nearest  neighbors  (1  As  and  1  Se). 

The  technique  was  also  used  to  investigate  impurity  charge 
shielding  arid  determine  atomic  arrangement  around  the  catalytic 
atoms  in  a  supported  catalyst^.  The  technique  is  finding  great 
interest  among  catalyst  chemists  at  a  time  when  progress  in  this  impor¬ 
tant  energy-related  field  is  in  the  national  interest. 

Accompanying  this  final  report  are  three  papers  (which  will  be 
submitted  for  publication)  describing  the  EXAFS  experimental^  and 
data  analysis  techniques^’  ^  These  reports  describe  in  detail  our 
practice  of  the  technique. 

For  the  future  of  the  technique,  a  proposal  was  made  to  NSF 
and  a  grant  awarded  to  the  University  of  Washington  (Stem,  Sayers, 

Lytle -consul  tan  L)  to  develop  an  EXAFS  experiment  at  the  St  an  fo  id  Synchro¬ 
tron  Radiation  Project.  The  x-ray  intensity  which  can  be  obtained  at 
this  source  is  at  least  10^  greater  than  that  of  ordinary  x-ray  tubes 
and  will  allcv,  data  of  higher  resolution  md  better  precision  in  a  much 
shorter  time. 

Conclusions 

In  competition  with  x-ray,  electron,  or  neutron  scattering 
the  EXAFS  technique  has  the  important  advantage  that  the  origin  of  the 
radial  structure  func'ion  (two  body  correlation  function)  is  fixed  at 
the  site  of  the  species-determined  x-ray  absorption  event  while  all 
other  techniques  average  over  all  the  different  kinds  of  atoms  in  the 
material.  The  results  obtained  on  amorphous  semiconductors  are 
important  and  unique  and  will  allcv  theorists  to  construct  more  accurate 
models  of  these  materials.  The  catalyst  work  is  only  the  top  o'  an 
iceberg.  This  was  the  first  example  of  a  technique  in  which  the 
environment  of  catalytic  atoms  in  general  could  be  determined. 

The  EXAFS  technique  has  been  shown  to  be  a  powerful  tool  for 
studying  the  arrangement  of  atoms  in  multi -component,  non- crys tall ine 
materials . 
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The  EXAFS  Technique:  I.  Experiment 
by 

-  *  * 

Farrel  W.  Lytle  and  Dale  E.  Sayers 

The  Boeing  Company,  Seattle,  Washington  98124 


Abs  tract 

Techniques  of  x-ray  spectroscopy  for  obtaining  extended  x-ray 

absorption  fine  structure  (EXAFS)  are  discussed.  The  essential 

instrumentation  is  an  automated  commercial  x-ray  diffractometer 

with  associated  power  supdy  and  counting  electronics  modified  to 

include  an  improved  x-ray  monochromator  alignment  fixture,  a  cryostat 

for  containing  the  absorber  sample  and  a  card  punch  for  data  readout. 

The  preparation  of  suitable  monochromator  crystals  is  discussed. 

A  preset  f’ount  mode  was  used  to  make  the  most  efficient  use  of 

experimental  time  and  the  technique  of  multi-pass  averaging  was  used 

to  obtain  good  counting  statistics;  i.e.  0.1%.  The  precision  and 

accuracy  of  the  wavelength/energ;  scale  was  determined  by  replicate 

analysis  and  comparison  to  emission  lines  from  the  x-ray  source  to 

be  —  1  leV.  When  these  lines  were  intense  the  nonlinear  response 

» 

of  the  x-ray  detecting  system  required  a  coincidence  correction 
to  be  evaluated  and  applied  to  the  measured  counting  rates.  EXAFS 
data  from  Cu  is  shown  as  an  example. 


*  Research  supported  in  part  by  the  Advanced  Research  Projects 
Agency  and  monitored  by  the  U.  S.  Army  Research  Office, 
Contract  No.  DAH  C04  71  C  0010. 
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Introduction 


This,  and  the  following  paper,  describe  the  experimental  details 
and  data  processing  techniques  which  we  have  developed  for  the  prac¬ 
tice  of  x-ray  absorption  spectroscopy  for  determination  of  the  struc¬ 
ture  of  non-crystalline  solids  (For  a  general  treatment  of  x-ray 
spectroscopy  techniques,  see  Sandstrom^.)  In  our  technique  the 
extended  x-ray  absorption  fine  structure  (dXAFS)  is  Fourier  analyzed 

to  produce  a  radial  distribution  function  centered  on  the  absorbing 

(2) 

atom  .  The  particular  utility  of  the  technique  is  the  ability  to 

explore  the  average  environment  of  each  atomic  species  in  multicom- 

(3) 

ponent  materials  such  aa  amorphous  semiconductors  or  the  highly 

(4) 

dispersed  materials  used  as  catalysts  .  The  basic  experimental 
determinant  of  the  technique  is  a  high  resolution,  precise  measure¬ 
ment  of  the  x-ray  absorption  coefficient  in  the  neighborhood  of 
characteristic  absorption  edges.  This  paper  describes  the  use  of 

a  commercial  diffractometer  (with  fairly  sample  modifications)  to 

« 

produce  an  x-ray  spectrometer  of  moderate  resolution;  the  diffrac¬ 
tion  crystals  and  preparation  of  the  diffractive  surfaces;  x-ray 
detectors  and  counting  techniques;  and  a  cryostat  and  sample  trans¬ 
port  mechanism. 

\ 

Summary  of  theory 

During  an  atomic  x-ray  absorption  event,  an  x-ray  with  energy 
ho  dissappears  and  a  photoelectron  of  kinetic  energy  E  is  emitted 
from  the  atom  conserving  energy  in  the  process 

E  -  ho  -  Ek  (1) 
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where  is  the  initial  binding  energy  of  the  electron.  The  absorp¬ 
tion  follows  the  law 

fix  -  lr.  Iq/I  -r  In  T/Tq  (2) 

where  p  is  the  linear  absorption  coefficient,  x  is  the  thickness, 

*6  is  the  incident  x-ray  intensity,  I  is  the  transmitted  intensity 
and  T  and  Tq  are  the  times  to  preset  count  with  the  absorber  in 
and  out  of  the  beam.  On  the  high  energy  side  of  the  absorption 

f 

edge  px  shows  fluctuations  (EXAFS)  which  are  related  to  the  atomic 
arrangement  immediately  surrounding  the  absorbing  atom.  We  have 
shown  that  the  EXAFS,  X(k)  can  be  accounted  for  by  a  scattering 
model ^ 

X(k)  *  1^2  E(r  )g  (2k)  exp(-2o  A2)  Sin[2kr  +  2r\  (k)]  (3) 

j  j  J  J  J  J 

where  k  is  the  electron  wave  vector,  is  the  number  of  atoms' in 

the  jth  coordination  sphere,  r^  is  the  radial  distance  to  the  jth 

atoms,  E(r^)  is  the  scattering  potential  encountered  by  the  electrons, 

gj (2k)  is  the  backward  scattering  amplitude,  the  exponential  is  a 

Debye-Waller  type  term  where  is  the  rms  fluctuations  of  the  atom 

about  r^  and  r)00  is  a  phase  shift.  The  form  of  this  equation  makes 

EXAFS  proportional  to  th.’  number  of  atoms  surrounding  the  absorbing 

atom  and  inversely  proportional  to  r^2.  A  decrease  in  temperature 

has  its  usual  effect  of  sharpening  -.tie  EXAFS.  Each  coordination 

sphere  contributes  a  sine-like  term  of  period  2kr  . 

Conceptually  EXAFS  may  be  visualized  as  a  kind  of  electron 
diffraction  where  the  source  of  -he  electrons  is  within  the  material 
originating  from  the  atoms  of  that  particular  atomic  species  involved 
in  the  abso-ption  event.  The  mechanism  of  the  absorption  fluctua¬ 
tions  is  as  »ollows:  during  the  dipole  transition  the  electron 
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oscillates  between  the  initial  and  final  states.  Since  the  final 
state  in  this  case  is  located  away  from  the  absorbing  atom  the  elec¬ 
tron  encounters  and  is  scattered  by  the  surrounding  potential  array. 
The  overlap  of  the  scattered  and  initial  wave  functions  in  the  dipole 
transition  matrix  produces  the  EXAFS  fluctuations. 
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Our  basic  x-ray  spectrometer  (as  shown  in  Fig  1)  is  a  hori¬ 
zontal  diffractometer  with  attached  tube  stand.  The  only  modifi¬ 
cation  is  an  improved  crystal  (monochromator)  support  for  optimum 
alignment  of  the  diffracting  crystal  (Fig.  lb).  This  spectrometer 
is  used  with  conventional  x-ray  diffraction  tubes  (Ag  cr  Mo  target) 
chosen  to  provide  a  high  intensity  continuum  in  the  region  of  the 
absorption  edge  to  be  measured.  The  problem  of  intense  spectral 
ines  which  sometimes  occur  if  discussed  in  the  section  on  counting 
techniques.  We  use  the  spectrometer  with  a  goniometer  radius  of 
18.75  cm  and  an  x-ray  tube  takeoff  angle  of  approximately  3°  (adjust¬ 
ed  for  maximum  diffracted  intensity).  The  first  slit  is  used  to 
limit  the  angular  divergence  of  the  x-ray  beam  (tube  spectrometer 
mode)  and  the  Bragg-Brentano  parafocusing  conditions  and  slit  posi¬ 
tioning  require  the  focusing  (exit)  slit  to  be  approximately  the 
same  size.  The  vertical  divergence  is  limited  to  ±  2°  by  a  Soller 
slit  located  before  the  exit  slit.  Thus,  the  resolution  function 
depends  primarily  on  the  size  of  the  entrance  slit.  (See  Klug  and 
Alexander^  for  a  complete  discussion  of  this  type  of  spectrometer.) 
Our  usual  method  of  operation  employs  0.05  mm  slits  (.025°  divergence) 
in  the  range  from  15°  to  45°  20  and  0.1  mm  (.05°  divergence)  slits 


for  angles  greater  than  45°  with  overlap  as  experimentally  required. 

Hence,  the  resolving  power  defined  as  0/AQ  is  approximately  1500. 

The  diffracted  intensity  from  the  continuum  operating  at  full  recom- 

3 

mended  tube  KVA  depends  on  wavelength  but  is  typically  1-10x10  cps. 

For  the  case  of  0.05  mm  slits  at  45°  20  and  a  LiF  crystal, 
the  energy  band  width  received  at  the  exit  slit  would  be  ~3.5eV. 
Since  the  intensity  distribution  of  the  radiation  filling  the  exit 
slit  has  the  usual  diffraction  profile,  the  Rayleigh  resolution 
criterion  suggests  that  the  spectrometer  should  be  advanced  in  angu¬ 
lar  increments  of  *5  the  angular  width  of  the  exit  slit. 

The  mode  of  ooeration  is  as  follows:  For  each  spectrometer 
position, T  and  Tq  are  measured  (preset  count)  and  stored  on  IBM 
cards,  the  spectrometer  is  advanced  to  the  next  20  position  and  the 
sequence  is  repeated.  A  separate  scaler  is  used  to  generate  a 
running  number  for  each  subsequent  pair  (T,  Tq)  .  Knowledge  of  the 
start  position  and  the  20  stepping  increm  nt  allows  calculation  of 
the  x-ray  wavelength  for  any  data  pair.  The  usual  experiment 
involves  500  -  1000  data  pairs.  The  spectrometer  stepping  motor 
*  and  absorber  changer  are  activated  and  synchronized  by  the  x-ray 
scaler  print-out  command. 

The  flanged  tube  stand  attached  to  the  spectrometer  allows  the 
x-ray  tube  to  be  transtated,  rotated,  and  inclined  with  respect  to 
the  spectrometer  circle.  The  alignment  procedure  consists  of 
locating  (by  means  of  these  adjustments)  the  most  intense  spot  on 
the  x-ray  tube  target  so  that  it  shines  through  the  zeroed  entrance 
and  exit  slits  as  measured  by  a  protected  x-ray  detector.  The 
crystal  monochromator  is  then  inserted  and  a  suitable  characteristic 


At  the  calculated  20 


line  chosen  for  final  crystal  adjustment, 
the  translation  and  rotation  ($)  adjustments  of  the  crystal  holder 
are  used  to  obtain  peak  diffracted  intensity  of  the  standard  line 
which  locates  the  diffracting  surfac.*  ot  the  crysta'  at  the  center 
of  the  spectrometer.  Further  refinement  of  the  alignment  is  no. 
necessary  if  the  procedure  described  in  the  section  on  precision 
and  accuracy  of  the  energy  scale  is  followed. 


Cryostat  and  Absorber  ^ 

EXAFS  has  been  shown  to  be  sensitive  to  temperature ^  \  The 

degree  of  the  effect  is  illustrated  in  Fig.  2b  and  2c  .which  were 

obtained  at  77°K  and  295°K,  respectively.  The  Debye-Waller  type 

exponential  occurring  in  the  theoretical  expression  for  EXAFS  can 

(8) 

be  related  to  the  usual  Debye  expression 


where  0D  is  the  Debye  characterisitic  temperature,  is  the  rms 
amplitude  of  vibration  of  the  atoms  in  the  jth  coordination  sphere, 

♦  K  is  Boltzmann  constant,  m  is  the  mass  of  the  jth  atoms ,  di  is 
B 

Planck's  ccnstant/2*.  $  076^  is  the  Debye  function,  and  T  is  the 
absolute  temperature.  Thus,  for  materials  with  a  low  value  of  0^. 
the  temperature  effect  will  be  largest,  however,  all  materials  will 
show  an  EXAFS  intensity  increase  (particularly,  far  from  the  absorp¬ 
tion  edge).  For  this  reasm  data  is  obtained  routinely  at  77  K 
(extra  sharpness  obtained  at  4  K  isn't  worth  the  effort).  It  is 


possible  to  use  the  temperature  effect  to  directly  evaluate  the  shell- 

(8) 

to-shell  coupling  coefficients 
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Figure  2.  Plots  of  In  T/  vs  the  energy,  E,  of  the  ejected 
photoe’.ectron  for  ..  2.5pm  thick  Cu  foil. 

Figure  2a.  One  pass  at  77°K,  105  preset  count.  The  features 

coincident  with  the  W  emission  lines  arc  due  to  the 
nonlinearity  of  the  x-ray  detection  system.  This 
effect  has  been  removed  on  the  other  graphs  using 
a  correction  technique  described  in  the  text. 

Figure  2b.  Ten  passes  at  77°K,  10^  preset  count  «.  xch  pass. 

Figure  2c.  Ten  passes  at  295°K,  105  preset  count  each  pass. 

Note  the  smearing  of  fine  structure  particularly 
at  high  energy. 


The  glass  cryostat  illustrated  in  Fig.  1  contains  5  liters 
of  liquid  nitrogen  which  will  last  for  3  days.  The  design  of  the 
sample  moving  mechanism  (Fig.  lh)  accomplishes  its  purpose  without 
movable  vacuum  seals  or  introduction  of  a  heat  leak  into  the  sample 
area.  If  a  stainless  steel  bellows  is  used  the  fatigue  life  is  10^ 
cycles  or  1-2  years  of  operation.  The  x-ray  transparent  windows  are 
1  mil  Mylar  expoxied  to  oval  openings  in  the  outer  shell. 
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The  absorber  samples  are  prepared  in  a  variety  of  ways;  malleable 
metals  are  rolled  (2-5jjm) ,  some  materials  may  be  evaporated  onto 
Mylar  or  thin  A1  foil,  soluble  materials  may  be  dissolved  and  then 
absorbed  and  dried  on  filter  paper,  many  materials  may  be  ground  as 
fine  as  possible,  mixed  with  a  suitable  plastic  glue  (Duco)  and 
solvent  (acetone)  then  cast  on  a  smooth  substrate.  When  dry,  the 
casts  are  sandwiched  between  \  mil  A1  for  support  and  thermal  contact 
and  attached  to  the  sample  holder.  The  optimum  absorber  thickness 
considering  contrast,  measurement  time,  and  primarily  the  coincidence 
counting  error  problem  is  attained  when  I/Iq~  1/3  on  the  high  absorp- 
ion  side  of  the  edge.  The  experiments  when  only  a  very  small  amount 
,  of  the  material  of  interest  is  present  must  use  thicknesses  such  that 
I/I  ~  1/10.  For  example  with  this  thickness  a  satisfactory  Ljjj 
absorption  pattern  of  0.3  wt  %  Pt  on  A^O^  was  ma<*e  *n  15  passes. 


Monochromator  Crystals 

The  usual  discussion  of  monochromators  for  x-ray  spectrometers 
emphasizes  resolution  with  narrow  crystal  rocking  curves  and  multi¬ 
crystal  spectrometers.  Our  requirement  stresses  high  x-ray  intensity 
at  moderate  resolution  for  good  statistical  accuracy  of  measured  EXAFS. 
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Although  better  resolution  would  probably  show  more  detail  there 

is  an  inherent  width  in  the  EXAFS  due  to  temperature  smearing  and 

(r  ) 

uncertainty  broadening  and  a  spectrometer  with  better  resolution 
than  the  natural  width  is  wasted  effort 

The  x-ray  intensity  is  limited  by  the  inherent  brilliance  of 
the  x-ray  tube  and  using  flat-crystal  optics  the  count  rate  is 
limited  by  the  efficiency  of  the  diffracting  crystal  and  the  width 

t 

of  the  energy  band  through  the  exit  slit.  The  diverging  x-rav 
beam  has  some  distribution  in  space  and  energy;  i.e.  photons/sec/ 
stearadian/eV.  The  monochromator  selects  with  some  diffraction 


efficiency  a  band  of  radiation  equal  to  the  angular  width  of  its 
rocking  curve  and  directs  it  toward  the  exit  slit.  The  optimum 
intensity  will  occur  when  the  divergence  of  the  incident  beam  as 
limited  by  the  first  slit  and  the  crystal  rocking  curve  are  appro¬ 
ximately  equal. 

The  efficiency  of  diffraction  and  the  rocking  curve  width  of 
diffraction  crystals  may  be  modified  by  appropriate  surface  treat¬ 
ment.  LiF  is  a  particularly  workable  crystal  in  this  respect.  In 
'Fig.  3  the  integrated  reflection  coefficient,  R,  for  non-polarized 
radiation  which  is  the  area  under  the  crystal  rocking  curve,  has 
been  calculated  as  a  function  of  wavelength  for  the  two  extremes  of 
crystal  perfection,  a  perfect  crystal  and  an  ideal  mosaic  crystal^) 
The  measured  values  of  selected  cleaved  and  treated  crystals  are 
also  shown.  Illustrated  in  Fig.  4  a  and  b  are  the  full  width  at 
half  maximum  (fwhm)  of  the  rocking  curve  and  the  diffraction  effi¬ 


ciency  (*diff racted/^incident^ 


at  the  peak  of  the  rocking  curve. 


The  measurements  were  made  on  a  double  crystal  spectrometer  in  the 


10 


X'  0.5 


>-'  0.4 

o 

z 

S 


!" 
<  •• 

cc 

u. 


•  CLEAVED 
A  ABRADED 
o  ETCHED 


MOSAIC 


WAVE  LENGTH.  A 


Figure  3. 


Integrated  reflection  coefficient  vs  wavelength  for 
LiF  as  cleaved,  roughened,  and  chemical  polished. 

The  measurements  were  made  on  a  double-crystal 
spectrometer  correcting  for  the  partial  polarization 
of  the  incident  x-ray  beam. 
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Figure  4b, 


o  .o 


The  intensity  ratio,  diffracted/incident,  at  the 
peak  of  the  diffraction  curve  for  LiF  as  a  function 
of  wavelength  and  surface  treatirent. 


WAVE  LENGTH.  A 


Figure  4a.  The  full  width  at  half  maximum  intensity  of  the 
rocking  curve  for  LiF  crystals  as  a  function  of 
wavelength  and  surface  treatment. 


(1,  -1)  mode  using  the  first  crystal  to  select  a  monochromatic 

beam  of  small  angular  divergence  (<  .01°).  The  test  crystal  was 

then  "rocked"  (rotated)  through  its  Bragg  angle  while  the  dif- 

(10)  , 

fracted  photons  were  measured.  This  is  a  direct  measure  of 
R  since 

R  *  £w_  M  f  JjLLLl  d0  (5) 

Iin.  Iin. 

9 

where  £  is  the  total  number  of  photons  counted  under  the  rocking 

curve,  w  is  the  angular  rate  of  crystal  rotation  (in  radians/  and 

I  and  I  refer  to  the  incident  and  diffracted  radiation, 

in.  dif. 

The  measurements  have  been  corrected  for  the  partial  polarization 
of  the  incident  test  beam. 

These  results  show  that  the  total  radiation  diffracted  by  LiF 

may  be  increased  by  factors  of  2-5  by  suitable  surface  treatment  to 

reduce  the  extinction  effects  which  occur  during  diffractions 

The  treatment  consists  of  vigorous  sanding  on  rough  paper  to  drive 

dislocations  into  the  crystal  followed  by  successively  finer  paper 

to  600  grit  to  provide  a  smooth  surface.  This  was  the  treatment 

'  o?  the  abraded  crystal  in  Figs.  3  and  4.  The  damaged  surface 

(12) 

layer  was  then  removed  by  an  etching  procedure  .  After  1  minute 
in  cone.  HF  a  chemical  polish  consisting  of  27.  NH^OH  in  H20  at  26  C 
with  vigorous  agitation  was  used  to  remove  the  surface  at  about 
lym/min.  The  crystal  was  repeatedly  checked  until  its  rocking  curve 
width  narrowed  to  that  of  the  slits.  Note  that  for  the  etched  crystal 
the  diffraction  efficiency  and  R  were  considerably  enhanced  over 
the  cleaved  crystal.  We  have  used  similar  techniques  with  Si  and 
quartz  crystals  with  some  success. 
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Counting  Technique 


The  precision  to  which  the  EXAFS  function  may  be  determined 
should  be  limited  by  the  usual  statistics,  og  3  where 
standard  deviation  and  N  -  number  of  photons  counted  for  each 
data  point.  This  is  usually  expressed  as  percent  counting  error, 

pee  «  (as/N)  x  100%;  e,g*  when  N  "  1C)6,  PCG  *  °,1S:*  W°uld 
appear  that  any  level  of  precision  could  be  obtained  by  making  N 

sufficiently  large;  but  in  practice  the  inherent  stability  (elec¬ 
trical  and  mechanical)  of  the  experiment  will  limit  the  precision. 
It  is  important  to  determine  this  value  for  a  particular  experiment 
otherwise  time  will  be  wasted  by  operating  at  a  larger  N  than 
warranted. 

When  the  inherent  experimental  precision  has  been  evaluated 
(as  described  below)  the  experiment  is  set  to  accumulate  the 
appropriate  number  of  counts  for  each  data  point.  If  better  pre¬ 
cision  is  desired  then  multi-pass  averaging  is  used  to  accomplish 


it. 


The  inherent  precision  of  cur  experiment  was  determined  as 

J 


follows:  For  a  range  of  N  from  10-10  a  subset  of  data  was  accumu¬ 
lated  for  each  value  of  N  and  3g  estimated  using^  standard  statistics 

(6)  ' 


Ji 

n  a. 


^  [K/t)-K/t)]2  ' 

_  m 


where  T  and  T  refer  to  the  time  to  accumulate  N  counts,  the  barred 
o 

value  is  the  mean  of  the  subset,  and  n  is  the  number  of  trials  in 
each  subset.  The  result  of  this  analysis  is  plotted  in  Fig.  5  where 
the  difference  between  the  measured  and  predicted  pee;  i.e. 

(measured)  -  N4]  100/N  is  plotted  vs  log  N.  Up  to  N  =  10 
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LOG  ACCUMULATED  COUNTS,  N 


Figure 


Figure  6. 


a. 


.  The  difference  between  the  predicted  and  measured 
percent  counting  error  vs  total  counts. 


lo/l  AND  CORRECTION  FACTOR 


Counting  rate  vs  absorption  ratio,  I  /I,  for  the 
detecting  systems  n,  b,  c  described  °in  Table  2. 

In  the  region  of  linear  response  the  ratio  is  constant 
then,  since  I  >  I,  1^  decreases  relative  to  I  and  the 
absorption  ratio  decreases.  Graph  d  is  a  plot  of  the 
correction  factor  (Short  's(l/')  method/  for  nyatemb 
which  can  be  multiplied  times  the  measured  count  rate 
to  correct  for  nonlinearity. 


the  data  shows  a  random  scatter  reflecting  the  accuracy  with  which 

a  was  measured  but  at  N  »  10^  there  is  a  sharp  break  in  the  ex^ri- 
s 

mental  cur^e  indicating  that  a  constant  level  of  precision  has  been 

attained  no  matter  h'  i  much  larger  N  was  made.  If  the  experiment 

had  been  completely  st.  .tistically  limited  the  pee  data  points  would 

cluster  around  the  zero  axis  for  all  N.  It  is  clear  from  this 

experiment  that  accumulating  more  than  10^  counts  (0.3%  pee)  is  a 

/ 

waste  of  time.  The  reason  for  this  fairly  high  pee  was  probably 
due  to  slight  mechanical  instability  in  the  crystat  support  which 
allowed  the  absorber  to  move  in  the  x-ray  beam  and  reflect  irregu¬ 
larities  in  thickness  in  the  measured  signal.  The  improved  pre¬ 
cision  apparent  in  the  EXAFS  of  Fig.  2b  and  2c  was  attained  by 

5 

averaging  ten  separate  runs  at  10  counts  for  each  data  point  for 

0 

a  final  precision  near  0.1%. 


Precision  and  Accuracy  of  the  Energy  Scale 

The  kinetic  energy  of  the  ejected  photoelectrons,  E,  must  be 
established  accurately  in  order  to  evaluate  the  natural  EXAFS 
variable  k  where 

k  - 

where  Aw  is  the  difference  between  the  work  function  and  the  aver¬ 
age  potential  seen  by  the  scattered  electron.  For  everv  experiment 
characteristic  and/or  impurity  lines  from  the  x-ray  tube  occur  and 
can  be  used  as  standard  reference  points  to  establish  the  energy 

(13) 

scale  at  the  accuracy  to  which  they  are  known  .  As  an  example, 
the  lines  indicated  in  Fig.  2  which  occurred  in  the  region  of  the 
experiment,  are  tabulated  in  Table  1.  Many  elements  are  present  as 


[ 


2m 

Xi* 


(E  -  Aw) 


(7) 


Table  1.  Precision  of  Energy  Scale 


Characteristic 

Line 

Wavelength. 

Measured  20 

Average  of  10  runs 

Calculated 

.  W  L  B9 

1.2C479+7 

34.838+3 

4.02459 

W  L  Bc 

1.21545+3 

35.162+2 

4.02396 

W  L  B2 

1.24460+3 

36.026+1 

4.02480 

w  l  e3 

1.26269+5 

36.564+1 

4.02523 

W  L  Bx 

1.281809+9 

37.142+2 

4.02478 

W  L  B6 

1.28989+7 

37.387+2 

4.02455 

W  L  64 

1.30162+5 

3_  733±2 

‘  4.0252/’ 

Cu  K 

1.392218+9 

40.473+2 

4.02496 

w  l  n 

] .42110+3 

41.350±2 

4.02502 

Avg. 4. 0247+3 


impurities  on  the  x-ray  tube  anode;  e.g.  W(always  present-eval- 

porated  from  the  filament),  Cu  (always  present-evaporated  from 

the  anode  block),  Pb,  Ni,  Ag,  and  Zn  are  usually  present,  plus 

the  characteristic  lines  from  the  primary  anode  element  in 

multiple  orders  of  diffraction.  The  presence  of  these  lines 

always  allows  a  standardization  to  be  made  as  veil  as  an  estimate 

of  the  precision  with  which  the  energy  variable  can  be  established. 

/ 

As  an  example  the  data  set  for  Fig.  2b  will  be  used.  Each 
of  ten  experimental  runs  produced  an  independent  value  of  2P  for 
each  standard  line.  The  ten  values  of  20  were  then  averaged  to 
obtain  one  value,  listed  in  Table  1  with  the  probable  error  deter¬ 
mined  from  the  replicatf  experiments.  Then  from  Bragg's  law  and 
the  tabulated  values(13),  2d  (twice  the  interplanar  spacing)  for 
the  monochromator  was  calculated  for  each  line  and  a  final  average 
obtained.  Following  this  procedure  for  every  set  of  data  removes 
the  requirement  for  corrections  involving  the  diffraction  process 
such  as  a  temperature  correction  for  the  monochromator  crystal, 
Lorentz-polar ization  correction,  refract<on  and  various  errors 
due  to  misalignment.  The  final  average  for  the  LiF(200)  monochro¬ 
mator,  2d  -  4.0247  ±  3,  is  a  standard  deviation  of  75  ppm  (typical 
of  our  experiments)  which  would  be  ±  0.7  eV  in  the  region  of  the  - 
Cu  K-edge  of  Fig.  2.  This  value  of  2d  is  used  in  a  computer  program 
to  calculate  (and  plot  as  in  Fig.  2)  the  value  of  E  associated 
with  each  data  pair  (Iq,  I)  since 

E  -  12398.04  eV.  (8) 

2d  Sin  0 

Note  that  the  error  reported  for  the  reference  lines  (~40  ppm) 
and  the  precision  to  which  they  can  be  measured  (~50  ppm)  are  comparable 


and  are  combined  in  the  final  determination  of  2d(  ~±  75  ppm). 

There  would  be  little  advantage  in  a  more  accurate  spectrometer 
unless  a  mere  accurate  reference  could  be  established. 

Coincidence  Losses 

In  Fig.  2a  artifacts  appear  in  the  absorption  curve  at  the 
position  of  the  strong  W  lines  W  L  6^  and  W  l  This  effect  is 

due  to  coincidence  losses  (nonlinearity)  occurring  in  the  x-ray 
detecting  system.  Although  predictable,  these  features  cause 
aerious  errors  in  EXAFS  data  analysis.  An  obvious  solution  is 
to  use  a  very  fast  system  with  a  large  dynamic  range;  however 
instances  will  occur  where  corrections  are  still  necessary. 

We  have  used  the  components  with  the  performance  capabilities 

described  in  Table  2. 

It  is  possible  to  obtain  good  data  with  any  of  these  systems 

(14) 

provided  corrections  are  made.  In  the  method  of  Short  the 
absorption  of  a  single  absorber  is  measured  as  a  function  of  counting 
rate.  When  the  system  becomes  nonlinear  the  measured  absorption 

4 

will  decrease  since  ux  *  In  Iq/I  and  Iq,  the  greater  countirg  rate, 
will  have  a  proportionally  greater  nonlinearity.  The  performance 
of  the  three  detecting  systems  is  shown  in  Fig.  6a,  b,  c.  Short  s 
technique  can  be  used  to  evaluate  a  correction  factor.  This  fretor 
is  plotted  in  Fig.  6d.  We  fitted  this  curve  with  an  equation  and 
during  computer  processing  all  data  was  corrected  to  remove  the 
effect  as  in  Fig.  2b  and  c.  The  correction  curve  must  be  evaluated 
to  an  accuracy  comparable  to  the  data,  i.e.  0.1 


IS 


PM  2106,  1.3ns  risetime. 


Discussion 


The  experimental  requirements  of  the  EXAFS  technique  are  within 
range  of  anyone  who  can  dedicate  a  commercial  x-ray  diffractometer, 
power  supply,  and  cotinting  rack  to  this  use.  The  necessity  for 
multi-pass  averaging  to  obtain  good  statistics  for  the  most  difficult 
£nd  most  interesting)  absorbers  requires  a  full-time  operation  with 
approximately  24  hours  for  each  pass  thro  ugh  ,  the  EXAFS  region.  Since 
the  experiment  is  easily  automated  this  is  net  a  handicap.  Indeed, 
our  experiments  are  turned  off  only  for  maintenance  every  few  months. 

The  additional  equipment  items  which  are  essential  are  1)  the 
improved  diffraction  crystal  holder.  (The  intensity  gained  by  opti¬ 
mizing  alignment  through  narrow  slits  can  be  a  factor  of  10.)  2)  The 

cryostat  and  sample  movement  device.  3)  A  card  punch  for  data  readout. 

Other  methods,  such  as  paper  tape  or  magnetic  tape,  lack  the 
advantage  of  immediate  visibility  and  easy  correction  of  occasional 
erroneous  data  points. 

After  the  data  is  acquired  and  the  multipass  average  is  obtained 
the  very  considerable  problem  of  data  analysis  remains  which  is  the 
'subject  of  the  following  paper. 
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Abs*  ract 

The  technique  of  computer  processing  of  extended  x-ray  absorption 
fine  structure  (EXAFS)  for  the  analysis  of  non- crystalline  solids  is 
discussed.  The  topics  include:  normalization  of  the  data  to  a  per  atom 
seal'  ,  Fourier  filtering  to  remove  the  sloping  background,  and  a  Fourier 
transform  of  the  data  into  a  rad.  structure  function.  Particular  factors 
peculiar  to  EXAFS  are  emphasized;  e.g.  upper  and  lower  limits  of  integration, 
the  effect  of  the  phase  shift  on  the  determin: tion  of  interatomic  distances, 
the  peak  shape  and  the  determination  of  coordination  numbers.  The  EXAFS 
spectrum  of  crystal  line  i*>  used  as  an  °-"unple. 
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Introduction 


The  previous  paper  described  the  experimental  requirements  for 
obtaining  EXAFS  data.1  This  paper  discusses  the  processing  and  analysis 
of  EXAFS  data  concentrating  on  the  important  factors  which  are  specifi¬ 
cally  related  to  EXAFS.  The  analysis  and  discussion  are  limited  to 
K-shell  absorption  spectra  only.  The  techniques  and  information  involved 
in  L-shell  spectra  will  be  published  elsewhere.2  The  general  subject  of 
finite  Fourier  transforms  has  been  much  discussed  by  others  involved  in 
x  ray  radial  distribution  analysis.3  This  paper  emphasizes  the  processing 
of  the  rav  data,  the  removal  of  the  EXAFS  from  the  average  absorption  back¬ 
ground  using  a  Fourier  filtering  technique,  and  the  normalization  of  the 
EXAFS  according  to  the  theoretical  expression  so  that  all  data  are  on  the 
same  per  atom  scale.  Fourier  transformation  using  the  point-scattering 
theory  results  in  a  radial  structure  function  centered  on  the  absorbing 
atom.  Tills  feature  permits  evaluation  of  unique  correlation  functions 
for  each  atomic  species  in  a  multi- component  material  and  is  the  major 
achievement  of  the  EXAFS  technique.  The  .articular  problems  of  changing 
from  an  energy  to  a  wave  vector  scale,  selecting  the  upper  and  lower  limits 
of  integration,  the  phase  Shift,  and  the  peak  shapes  are  discussed  in  detail. 


Initial  Data  Processing 


The  data  collected  from  our  h.XAFS  apparatus  consists  of  the  time 
to  collect  a  preset  number  of  counts  (usually  105)  with  the  absorber  in  and 
out  of  the  beam,  respectively  (T,Tq)  along  with  the  corresponding  step 
number  at  each  angular  position  of  the  spectrometer.  A  typical  data  deck 
consists  of  approximately  200-300  cards  with  4  sets  of  T,T  and  the 
running  number  per  card.  The  preliminary  data  processing  program  (1)  tab¬ 
ulates  the  initial  data  and  calculates  the  x-ray  energy,  the  kinetic 
energy  of  the  photoelectron,  and  the  total  absorption;  (2)  removes  the 
oscillatory  part  of  the  x-ray  absorption  from  the  smooth  monotonic  back¬ 
ground;  and  (3)  normalizes  the  data  to  a  per  atom  scale.  These  three 
program  functions  are  described  below. 


calculates  "h  K  "*  0,6  lnPU'  dMa  ““  P»*r-  routinely 

flln  V  ”  e“rgy  SC3le8  f“”  *“»«  data.  Speci- 

y’  ^  “Cal  1>  calculated  at  each  point  using 


yx  -  ln(T/T  ) 
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(l) 

*®  de3cribed  in  the  previous  paper.1  The  x-rav  0 

found  from  The  x  ray  energy  at  each  point  is 


hv  -  jj_398.04 

2d  Bl  i0 


(2) 


Where  (2d)  is  calculated  as  described  previously1  and  6  at  each  spectro- 

-er  set  1„8  ls  found  from  ^  ^  ^  ^  ’ 

"etlC  energy  °f  the  ejected  photoelectron  is  then 


E  -  hv-E, 


‘K 


(3) 


27  V  lnltlal  Mndla8  ener8y  °f  18  -«»•  which  -.-as  ejec-e 

Beard8  7 tSOrP“0n  °f  ^6  X'ra)r-  tabulated  binding  energies  of 

Bearden  and  Burr  were  used. 

The  next  step  in  the  processing  involves  removing  the  K  shell 

contribution  for  the  atomic  species  of  intprp^  f 

pecies  of  interest  from  all  other  absorption 

components  including  other  shells  in  the  sat*  atom,  absorption  by  other 
atoms  in  the  sample,  and  absorptim  by  the  specimen  support  which  is 

y  A1  foil.  In  the  energy  region  of  Interest  these  other  components 
y  fr‘>"  “y  characteristic  edges  and  may  be  adequately  des- 
cr  bed  by  a  Victoteen  formula7  of  the  form  „  .  ClW  where  A  is  the  x-ray 
wavelength  mad  C  and  0  are  fomd  by  fitting  this  formula  to  the  measured 
ab  orption  curve  over  an  energy  range  (typically  a  few  hundred  eV)  just 
below  the  K-edge  of  interest.  This  fitted  curve  is  then  extrapolated  to 
energies  above  the  absorption  edge  to  approximate  all  of  the  other  absorp¬ 
tion  components  and  is  subtracted  from  the  experimental  curve.  The 

mreT(k)  T  K  She11  “n"lbUtl0n  “  the  Including  the 

K 
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The  removal  of  the  oscillatory  part  of  the  x-ray  absorption 
coefficient  (EXAFS)  from  the  smooth.  »onotonlcally-de creasing  background 
absorption  Is  a  difficult  problem  since  neither  calculated  nor  empirical 
curves  exist  which  are  accurate  In  the  vicinity  of  the  absorption  edge. 

t  er  than  fitting  some  empirical  function  far  above  the  edge  and  extra¬ 
polating  back  to  the  edge  where  large  uncertainties  may  occur  w.  have 
Chosen  to  use  a  Fourier  filtering  technique.  This  technique  uses  a  fast 
curler  transform  algorithm  on  the  raw  data  (which  Is  equlspaced  In  9) 
and  takes  2048  data  points  and  analytes  them  In  terms  of  the  real  and 
imaginary  components  of  1024  frequencies.  The  lowest  frequency  has  a 
wavelength  of  twice  the  range  of  the  data  and  the  remaining  frequencies 
are  the  next  1023  harmonics  of  this  fundamental.  This  spectrum  has  3  com¬ 
pacts:  (1)  contributions  at  low  frequencies  due  to  the  monotonically 
decreasing  absorption  coefficient,  (2)  the  EXAFS  signal  at  Intermediate 
requencies ;  and  (3)  small  components  from  the  random  statistical  counting 
noise  uniformly  distributed  throughout  the  entire  frequency  range.  The 
lowest  EXAFS  frequencies  were  estimated  on  the  basis  of  the  typical  near 
neighbor  dlsr»nces  and  were  well  separated  from  the  low  freauencv  compon¬ 
ents  which  could  then  be  filtered  by  zeroing  the  first  few  low  frequencies 
and  re  transforming  the  spectrum.  This  Is  shown  in  Figure  1  where  70  of 
the  first  512  coefficients  of  a  typical  transform  are  sho™.  The  remaining 
coefficients  are  negligible.  The  arrow  Indicates  where  a  typical  cut 
would  be  made  to  filter  out  the  monotonlc  background.  This  transformed 
spectrum  contains  only  contributions  Iron  EXAFS  and  the  noise.  The  smooth 

absorption  background  may  be  found  by  subtracting  this  filtered  spectrum 
from  the  raw  data. 

The  entire  normalization  procedure  is  shown  in  Figure  2 
figure  2(a)  shows  the  Ge  K-edge  EXAFS  for  amorphous  GeSe  along  with  the 
extrapolated  Vlctoreen  curve  (dashed  line)  which  was  fitted  to  the  data 
n  the  euc-gy  region  below  the  edge.  In  Figure  2(b)  the  EXAFS  which  was 
obtained  from  the  data  in  Figure  2(a)  using  the  Fourier  filtering  technique 
la  shown  and  In  Figure  2(c)  the  smooth  background  curve  obtained  by  sub¬ 
tracting  tha  curve  In  Figure  2(b)  from  the  curve  In  Figure  2(a)  Is  also 
shown.  Figure  2(c)  also  compares  this  smooth  background  curve  u°  (k)  to 
both  a  Vlctoreen  formula  and  the  empirical  X2’73  curve  suggested^ 
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Figure  1.  Example  of  a  frequency  spectrum  of  the  raw  data  using  a 
fast  Fourier  transform.  The  separation  of  EXAFS  from 
the  low  frequency  slope  components  is  apparent.  The 
arrow  indicates  where  the  low  frequency  cutoff  is  made. 
Only  70  of  the  512  components  are  shown;  all  higher 
components  are  negligible. 
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Figure  2a.  The  measured  K-edge  of  Ge  in  amorphous  GeSe  where  px 
(absorption  coef.  times  thickness)  is  plotted  vs  the 
photoelectron  kinetic  energy,  E.  The  dashed  line  is  the 
Victoreen  formula  extrapolation  due  to  absorption  from 
other  than  the  Ge  K-shell. 

Figure  2b.  x  vs  E  normalized  and  extracted  from  Figure  2a  as  explained 
in  the  text. 

Figure  2c.  The  smooth  K-shell  absorption  background  y  (x)  vs  E  which 
was  used  to  normalize  x  according  to  equation  (A), 
y  (x)  is  obtained  from  Figure  2a  by  subtracting  the 
Victoreen  extrapolation  (dashed  line)  and  Fourier  filtering 
to  remove  the  EXAFS.  The  points  0  are  extrapolations  to 
the  edge  from  a  region  1000-1500  eV  above  the  edge  of 
both  a  Victoreen  form  and  the  empirical  formula  of  McMaster. 
Both  are  identical  and  approximately  5%  below  the  measured 


curve. 


McMaster  .  These  curves  were  normalized  to  the  data  in  the  region 
1000-1500  eV  above  the  edge  and  give  essentially  identical  results.  As 
can  be  seen  the  fit  is  good  to  500-1000  eV  but  is  lower  than  the  data 
curve  by  ^5%  near  the  edge,  which  would  give  noticeable  curvature  to  the 
EXAFS  (Figure  2(b)  if  they  had  been  used.  We  have  also  checked  the  sensi¬ 
tivity  of  this  technique  to  the  number  of  components  which  were  zeroed, 
and  have  found  no  significant  deviation  from  the  smooth  curve  as  long  as 
the  number  of  components  zeroed  was  less  than  the  expected  minimum  EXAFS 
frequency  (i.e.  if  no  more  than  the  first  10  frequencies  are  zeroed). 

The  final  step  in  the  initial  data  processing  involves  normal¬ 
izing  the  data  according  to  the  definition 

X(k)  -  £yR(k)/y°(k>J-l  (4) 

where  x(k)  is  the  EXAFS,  y^(k)  is  the  total  K-shell  contribution  to  the 

x-ray  absorption  coefficient  and  y°(k)  is  the  smooth  monotonically  decreasing 

*\ 

background  of  the  K  shell.  The  EXAFS  defined  in  this  way  normalizes  all 
data  to  the  same  per  atom  scale  which  will  allow  absolute  coordination 
numbers  to  be  determined  (as  described  below),  y  (k)  and  y  °(k)  are  deter- 
mined  as  described  above  and  the  normalized  EXAFS  x(k)  is  determined  using 
these  functions  and  equation  (4).  The  data  analysis  program  plots  yx  vs.  E 
and  x  vs*  E  and  punches  a  deck  of  x>E  pairs  which  is  used  as  the  input  for 
the  final  data  analysis. 


Fourier  Analysis  of  EXAFS 

Fourier  analysis  of  the  normalized  data  is  based  on  the  point- 
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scattering  theory  of  EXAFS  for  K-shell  absorption  ’  ’ 


X(k)  -  1/k  E 


jq  ■+ 

E  *2  E(rj)  gj(2k)exp(-20jk2)sin  ^2krj+2n(k)J  (5) 

J  rJ 


Here  g^(2k)  is  the  backward-scattering  amplitude  from  the  surrounding  atoms 
in  the  jth  shell  of  the  excited  p-state  electron  of  wave  number  k  per  unit 
electric  field.  The  direction  of  k  is  radially  out  from  the  excited  atom 
to  the  center  of  an  atom  in  the  jth  atomic  shell;  E(r^)  is  the  average 
electric  field  induced  at  rj(the  position  of  the  jth  shell  of  atoms)  by 


the  ionized  atom  at  r  ■  0.  N.  is  the  total  number  of  atoms  in  the  jth 
2  ^ 

shell  and  o^  is  a  measure  of  the  fluctuations  of  the  actual  positions  of 

the  atoms  in  the  jth  shell  about  their  average  position  r,.  The  exponen- 
2  7 

tial  containing  is  a  Debye -Waller-type  factor. 

Since  EXAFS  can  be  viewed  as  a  sum  of  damped  sine  waves  we  can 
formally  invert  the  data  using  Fourier  analysis.  In  most  monatomic 
materials  we  find  that  g(2k)  is  independent  of  r^  to  a  good  approximation 
so  that  equation  (5)  can  be  rearranged  and  transformed  with  respect  to 
sin  (2krf2n(k)  )  to  give 


<Kr) 


1 

\rn 


■^j-sin  (2kr+2n(k)  ) 


(6a) 


— L exp (-  (r-r  )2/2a2)  +  Ar 


(6b) 


where  equation  (6b)  is  obtained  from  equation  (6a)  using  equation  (5),  and 
$  (r)  is  a  radial  structure  function  which  given  iha  Icccticr.  of  ctcms 
around  the  absorbing  atom  in  terms  of  a  series  of  gauss ian-sh aped  peaks 
whose  amplitudes  are  proportional  to  N^.  Ar  contains  a  symmetrical 
exponential  tarm  which  is  not  used  find  other  terms  much  smaller  than  <Kr). 


While  the  integration  of  equation  (6)  is  illustrative  to  show 
the  kinds  of  information  available  from  EXAFS  it  ignores  several  factors 
which  must  be  accounted  for  in  the  actual  data  analysis  including  the  finite 
upper  and  lower  limits  of  integration,  the  fact  that  g(2k)  and  n  (k)  are  not, 
in  general,  known  and  that  the  definition  of  k  is  uncertain.  Specifically 
the  integral  that  is  actually  calculated  is 


$(r) 


vrr  / 

2tt  k 


max 


min 


x(k)  exp  (2ikr)dk 


(7) 


The  problems  involved  in  equation  (7)  fall  in  2  categories 
(evaluation  and  interpretation)  and  will  be  discussed  separately. 


A)  Evaluation  -  The  photoelectron  wavevectcr  is  evaluated  from  the 
definition 


—  "  E  "  (8) 

where  E  is  defined  by  equation  (3)  and  is  referenced  to  the  Fermi  energy, 

Uo  is  the  average  inner  potential  measured  from  the  vacuum  level  and  W  is 
the  work  function.  For  metals  Iq-W  is  on  the  order  of  the  Fermi  energy 
measured  from  the  bottom  of  the  conduction  band.  In  semiconductors  and 
insulators  the  zero  of  kinetic  energy  is  not  as  easy  to  define  but  may  be 
approximated  by  averaging  the  pseudo-potential  over  the  energy  range  of 
interest.  In  either  case  errors,  AE,  of  a  few  eV  may  arise  in  the  choice 
of  Uq-W  used.  These  errors  do  not  significantly  affect  the  data  since 
even  at  the  energy  corresponding  to  the  lower  limit  of  integration  (typi¬ 
cally  El  ^50eV),  AE/El  is  negligible.  We  have  verified  this  on  several 
spectra  by  varying  Uq-W  over  a  suitable  range  and  found  no  significant 
effects  on  the  transforms. 

Another  problem  in  evaluation  is  the  cnoice  of  upper  and  lower 
limits  of  integration  in  equation  (7).  The  upper  limit  is  determined 
primarily  by  the  experimental  range  of  energies  which  have  been  measured, 
typically,_1000-1500  eV  above  the  edge  (which  corresponds  to  a  k  range  up 
to  15-20  A  )  unless  other  absorption  edges  appear  in  the  data.  It  is 
usually  desirable  to  have  k  as  large  as  possible  to  reduce  the  termination 
error  found  in  Fourier  transforms  which  are  truncated3  at  finite  k,  although 
these  effects  are  not  the  primary  cause  of  deviations  from  the  gaussian- 
peak  shape  predicted  by  equation  (6b)  as  will  be  discussed  below.  The 
lower  limit  of  integration  is  generally  determined  by  the  requirement  that 
it  is  far  enough  above  the  edge  so  that  the  assumptions  of  the  point 
scattering  theory  are  valid  and  any  band  structure  effects  on  the  density 
of  states  have  disappeared.  There  is  no  precise  point  at  which  this 
occurs  but  is  typically  .30-50  eV  and  may  be  verified  by  choosing  several 

lower  limits  and  selecting  the  one  which  produces  the  cleanest  transform 
(i.e.,  low  noise  at  r<l£). 


The  actual  transform  may  be  done  either  by  direct  integration 
or  by  using  a  Fast-Fourier  algorithm.  In  order  to  process  large  amounts 
of  data  with  the  capability  to  make  real  time  decisions  we  have  developed 
a  program  using  a  display  scope  with  interactive  capability  .  The  out¬ 
puts  are  the  real  (cosine),  imaginary  (sine),  and  magnitude  of  the  trans¬ 
form.  In  Figure  3  the  EXAFS  for  crystalline  Ge  over  the  range  to  be  trans 

formed  is  shown.  In  Figure  4  the  transform  of  the  data  from  Figure  3  is 

o-l 

shown  over  the  range  of  k  from  4.1  to  15. SA  ,  where  4(a)  is  the  real  part 
4(b)  the  imaginary  part,  and  4(c)  is  the  nagnitude  of  the  transform.  As 
can  be  seen  the  imaginary  and  magnitude  ha  re  several  peaks  which  contain 
structural  information.  The  factors  which  «  ffect  the  interpretation  of 
these  transforms  and  the  information  contained  in  the  transform  is 
discussed  in  the  next  section. 

B)  Interpretation  -  There  are  several  qualitative  features  which  help 
understand  the  data.  First,  as  can  be  seen  in  Figure  4,  there  is  very 
little  amplitude  in  the  region  below  1A.  This  condition  can  usually  be 
achieved  by  proper  selection  of  the  lower  limit  of  integration  and  indi¬ 
cates  that  there  are  no  large  normalization  or  background  subtraction 

errors  which  have  been  Vnown  to  affect  the  transforms  strcngly  in  ..his 
11  12 

region  '  .  Because  of  the  uncertainties  which  have  been  discussed  in 

removing  the  EXAFS  from  the  background  absorption  it  is  likely  that  some 
error  has  been  introduced  in  the  transform;  however,  our  experience  is 
that  these  errors  can  be  made  small. 

The  relationship  between  the  real  part,  imaginary  part  and 
magnitude  of  the  transform  is  apparent  in  the  first  large  peak  in  the 
data  (between  1-3A) .  If  the  theoretical  expression  from  equation  5  is 
inserted  into  equation  7  and  it  is  assumed  that  n”0  and  the  k  dependence 
of  the  factor  g  (2k)/k  is  neglected  the  result  would  be  a  symmetric 
gauss ian-shaped  imaginary  part  and  an  antisymmetric  real  part  centered 
on  the  shell  distance,  r  .  This  is  what  is  observed  in  Ge  for  the  first 
shell.  The  peak  shape  of  the  imaginary  part  has  deep  negative  side 
lobes,  which  will  be  discussed  later,  but  this  does  not  change  the  basic 
relationship  between  the  two  parts  of  the  transform.  In  fact,  the  two 
parts  of  the  transform  contain  complementary  information  so  that  only 


y(k)  vs  k  (A  )  for  crystalline  Ge.  Data  is  shown  on  .y  in 
the  region  4.1  -  15.8  X"1  which  is  Fourier  transformed  in 
Figure  4. 
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Figure  4a, 
4b,  4c. 


The  real  part,  imaginary  part,  and  magnitude,  respectively, 
of  the  Fourier  transform  of  the  EXAFS  in  Figure  3  plotted 
4>(r)  vs  r(&) .  The  numbers  over  the  peaks  locate  the  first 
6  coo  rd  in  at  i  o.i  shells.  The  $(r)  scale  in  4c  is  twice  that 
of  4a  and  4b. 


one  part  need  be  studied.  The  situation  is  somewhat  complicated  by 
the  finite  constant  phase  shift.  By  using  simple  trigonometric  iden¬ 
tities  it  can  be  shown  that  the  real  and  imaginary  parts  now  are 
admixtures  of  the  rpO  transforms  so  that  neither  part  is  completely 
symmetric  or  antisymmetric  unless  the  phase  is  a  multiple  of  tt/2. 
Location  of  the  peak  position  may  still  be  made  using  the  magnitude  of 
the  transform  which  is  independent  of  any  constant  phase  shift  and  will 
be  peaked  at  r'srj*  In  the  case  of  Ge  it  appears  that  the  constant  part 
of  the  phase  is  approximately  zero  since  the  real  aid  imaginary  parts 
of  the  Fourier  transform  have  the  shapes  described  above.  Note  that 
the  magnitude  (Figure  4c)  is  peaked  at  the  same  place  as  the  imaginary 
(Figure  4-)  part  but  is  smoother  and  broader  which  in  some  ways  makes 
it  easier  to  identify  shell  locations  since  it  reduces  the  effects  of 
the  side  lobes. 

To  locate  the  peak  position  accurately  it  is  necessary  to 
understand  the  phase  shift  more  fully.  Estimates  of  the  phase  shifts 
using  the  WKB  ^  method  have  shown  that  for  most  elements  the  phase 

shift  is  approximately  linear  in  k  over  the  EXAF3  energy  range  with  a 
negative  slope,  i.e. 


n(k)  -  -ak+0 


(9) 


with  a,$  constant.  This  linearity  has  been  verified  in  the  case  of  Cu 

using  a  more  exact  calculation  of  phase  shifts  using  a  Hartree-Fock 

potential.1"*  Substituting  equation  (9)  into  equation  (5)  gives  the 

argument  for  the  sine  wave  of  (2kr-2ak+2B ) .  The  effect  of  the  constant 

phase  2$  has  been  discussed  above;  however,  as  first  pointed  out  by  Mott1^, 

the  linear  k  dependent  term  shifts  the  "frequency"  of  the  sine  wave  from 

ij  to  r^-a.  In  the  Fourier  transform  this  has  the  effect  of  shifting 

all  peaks  toward  the  origin  by  an  amount  a.  In  crystalline  materials 

where  the  distances,  r  ,  are  known  the  transforms  may  be  used  to  evaluate 

a.  This  is  shewn  in  Table  I  where  the  data  for  the  transforms  shown  in 

Figure  4  is  tabulated.  Column  five  of  the  table  shews  that  for  the  first 

o 

3  shells  where  the  data  is  most  reliable  a  value  of  a  »  .28A  may  be  used. 


Table  I.  EXAFS  Structural  Data  For  Ge 
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tjCA) 

*  0 
rj(A> 

0 

a  (A) 

y*i 

1 

4 

2.45 

2.17 

.28 

i 

2 

12 

4.00 

3.72 

.28 

.45 

3 

12 

4.68 

4.41 

.27 

.30 

4 

6 

5.65 

5.38 

.27 

.04 

5 

12 

6.16 

5.83 

.33 

.06 

6 

24 

6.92 

6.60 

.32 

.13 
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known  jth  shell  distance  for  dia.  cub.  Ge  (d»5.65A) 
measured  peak  position  from  Figure  4b 


“  rj"rJ 


Aj/A^  -  amplitude  of  the  imaginary  transform  at  r^  normalized 
to  the  first  shell  amplitude 


This  data  also  shows  at  least  6  shells  in  Ge  may  be  located  although 

o  o  t 

there  are  other  features  in  the  data  at  ^5A  and  ^6.2A  which  can  t  be 
explained  by  the  simple  point  scattering  model.  The  value  of  a  which^ 
can  be  determined  from  known  structures  is  typically  accurate  to  ±.02A 
which  gives  a  total  accuracy  of  2-3%  in  the  distance  measurements.  To 
determine  distances  in  materials  where  the  crystal  structure  is  not 
known  (such  as  amorphous  materials  or  biological  systems)  assumptions 
must  be  made  about  the  phase  shift.  Generally,  the  phase  shift  is  domi¬ 
nated  by  the  potential  well  created  by  the  inner  vacancy  in  absorbing 
atom  itself  and  not  by  its  surroundings  so  that  all  Ge  atoms  should  have 
approximately  the  same  phase  shift  regardless  of  the  compound  it  is  in. 

This  is  verified  by  the  data  in  Table  II  where  the  k  dependence  of  the 
phase  shifts  in  a  variety  of  materials  whose  structures  are  known  are 
listed.  The  Ge  phase  shifts  in  the  Ge  compounds  agree  very  well  with  the 
value  of  a  *  .28  determined  from  crystalline  Ge .  The  values  of  a  for 
several  other  elements  in  the  third  row  of  the  periodic  table  are  listed 
and  show  that  a  is  between  .25  and  .41  in  this  region  of  the  periodic 
tabic.  To  minimize  the  uncertainties  in  assigning  a  for  unknown  structures 
we  have  used  standard  compounds  whose  composition  and  local  structure  are 
known  and  whose  local  environment  is  expected  to  be  similar  to  the  local 
structure  being  investigated.  As  a  general  rule  we  have  found  that  com¬ 
pounds  ancng  atoms  in  the  same  row  of  the  periodic  table  (especially  if 
the  valence  of  the  absorbing  atom  is  small)  tend  to  have  phase  shifts  similar 
to  the  pure  element  whose  edge  is  being  studied.  We  have  also  found  that 
in  some  systems,  particularly  certain  oxides  with  ionio  bond,  there  are 
significant  deviations  from  the  elemental  phase  shifts.  Further  experi¬ 
mental  and  theoretical  studies  are  underway  in  order  to  quantitatively 
understand  the  observed  phase  shifts  and  their  origins  and  make  reliable 
structural  measurements  in  a  wider  variety  of  systems. 

As  mentioned  above  the  shapes  of  the  peak  differ  significantly 
from  the  gaussian  shape  predicted  by  equation  (6b)  even  if  termination 
effects  are  taken  into  account.  This  is  because  in  the  actual  transform 
(equation  (7))  the  k  dependence  of  the  factor  g(2k)Ac  has  not  been  removed 

from  the  data  so  that  the  experimental  transform  is  a  convolution  of  the 


Table  II.  Phase  Shift  k  Dependence  (g  from  equation  9) 
for  Several  Elements  and  Compounds 


Material 

a 

Cr 

.36 

Fe 

.26 

Ni 

.25 

Cu 

• 

CO 

00 

Zn 

.41 

As 

.27 

Mat  e  ri  al 

a 

Ge 

CO 

CM 

• 

Se 

.26 

(Ge)02 

.27 

(Ge)Se2 

.29 

(Ge)Se 

.29 

result  of  equation  (6b)  with  the  Fourier  transform  of  g(2k)/k.  The 
general  shape  of  this  function  is  the  envelope  of  the  EXAFS  shown  in 
Figure  3  which,  for  many  systems,  increases  to  a  peak  about  k  =  6-7  A*" 
and  decreases  at  higher  k.  This  envelope  also  includes  temperature 
effects  but  these  decrease  monotonically  with  k  and  do  not  give  the  peaking 
which  is  observed.  The  transforms  of  a  variety  of  similar  envelope 
functions  have  been  studied  by  Waser  and  Schomaker^  who  show  that  the 
results  are  transforms  with  large  negative  side  lobes  as  observed  in 
Figure  4.  At  present  the  functional  form  of  g(2k)  is  not  known  well 
enough  to  separate  its  effect  from  the  data.  This  complicates  the 
analysis  particularly  at  large  r  where  the  shells  are  close  enough 
together  so  that  the  side  lobes  interfere  with  adjacent  peaks  making 
peak  identification  uncertain.  The  peak  shape  also  makes  interpretation 
of  the  amplitude  difficult.  From  equation  (6b)  the  amplitude  of  each 
peak  is  proportional  to  E(  r j  ^ j  r j *  by  a  standard  integration  the 

area  under  the  peak  is  proportional  to  N^r^/r^.  However,  the  area  is 
difficult  to  measure  accurately  so  the  amplitude  is  usually  used  to 
obtain  any  information  about  number  of  atoms  per  shell.  Since  r^  is  deter¬ 
mined  from  the  peak  location,  the  amplitude  contains  information  propor¬ 
tional  to  NjE(r^)/Oj.  When  the  data  is  normalized  according  to  equation  4, 
all  data  are  on  the  same  scale  and  could  be  compared  directly  to  deter¬ 
mine  the  coordination  numbers,  NJf  if  E(rj)/oj  were  known.  Normally,  this 
is  not  the  case  so  reliable  coordination  numbers  are  obtained  by  empirical 
comparisons.  In  some  crystalline  systems  o^  may  be  calculated  and  the 
measured  amplitude  used  with  the  known  crystal  structure  to  determine 
E(rj);  however,  it  is  usually  not  possible  to  measure  enough  amplitudes 
reliably  to  differentiate  between  possible  forms  for  E(rj)  (e.g.  power  law 
vs.  exponentail) .  It  is  also  possible  to  measure  the  temperature  dependence 
of  the  amplitudes  to  determine  o^  experimentally.  In  multi-component 
materials  where  only  one  absorption  edge  can  be  measured, ^  we  have  assumed 
that  E(r^)  is  primarily  a  function  of  local  environment  and  used  standard 
compounds  whose  environment  and  coordination  are  expected  to  be  similar  to 
the  unknown  system  to  determine  N^.  We  estimate  that  this  procedure  can 
give  coordination  numbers  which  are  good  to  >v 207.  for  the  first  few  shells. 


L 


- - ~ 
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Recently,  we  have  extended  this  procedure  to  making  absolute  measure- 

ments  of  the  first  neighbor  coordination  of  both  species  in  binary  com-  ^ 

pounds  where  each  atom’s  EXAFS  can  be  measured.18  It  was  possible  to 
determine  the  number  of  each  kind  of  first  neighbor  atom  around  both 
species.  This  technique  was  limited  to  systems  for  which  good  standards 
could  be  used  for  a  reference,  but  did  give  reliable  results  for  amorphous 
GeSe2  and  As^Se^  which  are  similar  to  their  crystalline  polymorphs.  To 
obtain  more  reliable  coordination  numbers  will  require  a  better  under¬ 
standing  of  the  meachanisms  which  determine  E(r)  and  g(2k),  so  that  these 
factors  can  be  determined  reliably  without  the  use  of  standard  compounds. 


Discussion 


The  Fourier  transform  of  EXAFS  data  can  be  used  to  determine 
the  local  atomic  environment  of  each  kind  of  atom  in  multi- component 
non-crystalline  materials.  It  is  a  general  and  powerful  tool  for  study¬ 
ing  any  physical  state  of  matter,  solid,  liquid,  or  gas.  It  has  been 
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applied  to  amorphous  semiconductors,  ’  ’  supported  catalysts  and 
biological  materials  (yet  to  be  published) . 

The  data  processing  is  concerned  with  extraction  of  the  EXAFS 

from  data  containing  the  effects  of  other  physical  processes  in  a  form 

suitable  for  Fourier  analysis,  which  is  a  difficulty  common  to  other 

3 

scattering  methods.  However,  there  are  particular  factors  peculiar  to 

EXAFS  such  as  the  phase  shift,  definition  of  the  wave  vector  scale,  and 

an  unknown  scattering  factor  which  are  the  primary  concerns  of  this 

paper.  We  have  addressed  these  problems  in  a  completely  pragmatic  way 

by  using  well-matched  standard  materials  to  determine  these  parameters 

and  then  applying  them  to  the  unkncwn  materials.  As  such,  the  technique 

is  useful  for  studying  the  atomic  arrangement  in  non-crystalline 

materials  but  suffers  a  lack  of  completeness.  As  more  materials  are 

studied  and  the  parameters  become  more  accurately  defined,  we  trust  that 

theory  will  catch  up.  In  addition  to  the  structural  information 

emphasized  here,  the  x-ray  absorption  edges  contain  a  wealth  of  physics 

which  is  becoming  apparent  as  we  have  tried  to  understand  the  EXAFS 

technique;  e.g.  some  L-spectra  show  giant  resonances  at  the  edge  with 

2 

s-d  mixing  in  the  EXAFS  and  the  role  of  impurity- charge  shielding  in 

9 

EXAFS  has  been  discussed. 
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ABSTRACT 


A  computer  program  is  described  which  allows  real-time 
data  analysis  decisions  by  means  of  a  display  scope  for 
Fourier  transforms,  data  manipulation,  and  inverse  trans¬ 
forms  of  x-ray  data.  The  example  given  analyzes  the  Ge 
K-edge  extended  x-ray  absorption  fine  structure;  however, 
any  x-ray,  electron,  or  neutron  scattering  data  could  be 
analyzed. 


INTRODUCTION 


Fourier  transforms  of  experimental  data  are  widely  used  in  many  areas 
of  physics  including  x-ray  and  electron  scattering  and  extended  x-ray 
absorption  fine  structure  (EXAFS)  experiments  as  well  as  many  other 
applications.  Often  these  transforms  involve  large  amounts  of  data 
and  are  very  time-consuming;  however,  the  development  of  the  Fast 
Fourier  Transform  (FFT)  algorithm  (Cooley  &  Tukey,  1965)  resulted  in 
substantial  time  savings.  Data  analysis  may  still  be  time-consuming, 
particularly  if  real-time  decisions  must  be  made  regarding  optimum 
termination  points,  multiplication  by  termination  functions,  inverse 
transforms,  and  display  of  final  data,  which  traditionally  means  several 
individual  computer  runs  with  resultant  turnaround  time. 

These  shortcomings  have  been  overcome  with  the  development  of  a  data 
analysis  proqram  which  uses  a  vector  display  scope  on  line  with  a  large 
digital  computer.  This  interactive  computer  program  allows  rapid  initial 
processing,  real-time  decisions  regarding  the  data  with  a  light  pen  or 
typewriter,  and  FFTs  all  of  which  require  only  a  few  minutes  per  data  case. 
In  the  following  sections  the  hardware,  software  and  a  sample  data  analysis 
of  EXAFS  will  be  used  to  demonstrate  the  program's  capabilities. 

The  program  has  been  developed  to  analyze  digital  EXAFS  data  (see  Sayers, 
Lytle,  &  Stern,  1971  for  examples),  but  can  be  used  for  any  other  kind  of 
Fourier  analysis  problem.  Data  in  graphical  form  may  also  be  analyzed  by 
digitizing  first.  Most  scattering  experiments  (x-ray,  electron,  neutron 
or  EXAFS)  are  conducted  so  that  data  are  accumulated  in  equi-spaced  angular 
increments  which  are  not  equally  spaced  in  the  transform  variable.  This 
interpolation  is  accomplished  as  a  first  step  in  this  program,  hence  any 
data  equi-spaced  or  not  may  be  analyzed. 


1.  COMPUTER  HARDWARE  ' 


The  hardware  for  this  program  is  a  System/360  model  65  IBM  computer, 
attached  directly  to  an  IBM  2250  vector  display  scope  with  a  32-key 
function  keyboard,  a  standard  typewriter  type  entry  keyboard,  and  a 
light  detecting  pen.  The  scope  is  similar  in  appearance  to  a  TV 
screen.  The  light  detecting  pen  can  be  touched  directly  on  the 
screen  yielding  the  precise  coordinates  of  contact,  which  can  be 
used  both  for  terminating  the  data  and  issuing  program  instructions. 
Internal  core  storage  is  limited  to  280K  bytes  at  our  facility  because 
of  the  mul ti -processing  requirements  of  the  computer.  This  problem 
has  been  reduced  by  development  of  an  FFT  algorithm  (Singleton,  1967) 
which  uses  a  relatively  small  amount  of  in-core  storage  and  uses  10 
auxiliary  mass  storage  units  for  data  storage.  The  input  data  is 
entered  on  cards  with  the  program  residing  on  permanent  disc  storage. 
Printed  output  goes  on  a  standard,  high-speed  printer.  Graphical 
output  (hard  copies)  of  any  screen  display  may  be  obtained  on  either 
SC4020  or  CALCOMP  plotting  hardware  at  the  user's  option.  This  output 
Is  written  on  two  7-track  tapes  to  be  processed  offline.  The  entiie 
program  requires  280K  bytes  of  internal  mass  storage  and  operates  well 
in  the  multi-variable  task  processing  environment. 

2.  COMPUTER  SOFTWARE 

The  interactive  computer  graphics  program  consists  of  one  long  driver 
program  that  calls  numerous  subroutines  from  the  IBM  2250  graphics 
subroutine  package  (IBM,  1968)  and  the  FFT  subroutine.  All  mathematical 
manipulation  with  the  exception  of  the  transform  is  performed  by  the 
main  program.  The  IBM-supplied  software  facilitates  communication  with 
the  main  program  in  controlling  light-pen  and  function  key  interrupts, 
display  regeneration,  and  typed  entries. 

At  our  facility  core  storage  for  scope  jobs  is  limited  to  280K  bytes  of 
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core  which  'represents  a  problem  since  Fourier  transforms  of  2  complex 


(Continued) 

data  pairs  must  be  performed.  The  graphics  subroutine  package  takes 
about  170K  bytes  of  storage,  leaving  only  110K  bytes  for  the  main 
program  and  the  FFT  algorithm.  Generally,  FFT  algorithms  (Bingham 
4  Godfrey,  1967)  need  at  least  2.5  times  the  sue  of  the  trans  orm 
for  internal  storage  in  order  to  perform  an  in-core  transfora.  Sin 
this  amount  of  storage  was  simply  not  available  for  * 

graphics  mode,  an  algorithm  was  chosen  that  uses  auxiliary  mass  st  rage 
Disc  Files)  (Singleton,  1967).  This  FFT  algorithm  occupies  only 
bytes  of  storage  plus  20K  bytes  for  high-speed  buffers  to  per  form  a 

double  precision  transform.  All  of  the  routines  used  inclu  *  ' 

AnA  in  pfiRTRAN  with  the  exception  of  certain 
FFT  subroutine  were  coded  in  FUKIKAIm  wun  u.c  p 

GSP  routines. 


ANAIYSIS  PROCEDURE 

order  to  demonstrate  the  capabilities  of  this  program,  a  sample  data 

examination  is  presented  using  the  K-edge  EXAFS  from  crysta  ine 

/r  \  rvAcc  is  the  series  of  oscillations  appearing  on  the 
germanium  (Ge).  EXAFS  is  the  ser  (Sayers, 

high  energy  side  of  x-ray  absorption  edges.  The  theory 

et  al.,  1971)  is  given  by 


„(■„).!  Nd  g(k)  exp(-c^k2)  Sin(2krj  *  2n) 

j=l  r? 


hniic  nf  N  atoms  surrounding  the  absorbing 
wherG  *“  iS  ""  ‘heq(Sk  1  In  dependence  of  the  scattering 

atom  at  a  distance  r.,  g(  oebye-Waller  type  temperature 

of  the  photo-electron,  the  exponential  is  a  Debye  uai  yu 

-  -  *  *  ^  *-vmd ’t— r  nr  in  r  nr 

areas  are  proportional  to  This  technique  has  bftn?appl i  ^ 

studies  of  amorphous  semiconductors  (Sayers,  e  a  & 

(Lytle,  et  al.,  1974),  and  charge  shielding  in  me  a 

1973). 
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1.  DISPLAY  CHI  VS  K  FOR  THE  DATA  JUST  READ  IN. 

2.  CHOOSE  LIMITS  IN  K  SPACE  FOR  FOURIER  TRANSFORM. 

3.  DISPLAY  CHIBAR  VS  KBAR  FOR  THE  RANGE  TO  BE  TRANSFORMED 
AFTER  THE  DATA  HAS  BEEN  EOUISPACED  IN  K. 

4.  MULTIPLY  BY  CONVERGENCE  FACTOR,  A-EXP(-B-k2) . 

5.  PERFORM  THE  FOURIER  TRANSFORM. 

6.  DISPLAY  THE  TRANSFORM: 

REAL 

IMAGINARY 

MAGNITUDE 

7.  MULTIPLY  THE  TRANSFORM  BY  Cr\ 

8.  SELECT  THE  LIMITS  FOR  THE  INVERSE  TRANSFORM. 

9.  PERFORM  THE  INVERSE  TRANSFORM. 

10.  DISPLAY  THE  INVERSE  TRANSFORM. 

11.  READ  THE  NEXT  DATA  CASE. 

FIGURE  1  -  "Menu"  of  Program  Options 


(Continued) 


One  or  more  decks  of  Chi  versus  E  (E  is  the  energy  above  the  x-ray 
absorption  edge)  arc-  read  into  the  computer.  The  program  displays 
the  "menu"  (shown  in  Figure  1)  on  t,  e  scope  listing  the  analysis 
options  available.  Normally  analysis  proceeds  step-by-step  return¬ 
ing  to  the  menu  after  each  step;  however,  if  desired,  any  step  may 
be  repeated.  The  options  are  shown  in  the  flow  chart  in  Appendix  I 
and  are  selected  by  touching  the  light-pen  on  the  line  of  the  menu 
desired. 

Normally,  the  first  operation  is  to  display  the  input  data  with  a 
change  of  independent  variable  from  energy  (E)  to  the  electron  wave 
vector  k  which  is  the  appropriate  Fourier  variable  and  is  related 
to  E  by 


where  m  is  the  electron  mass,  Ai  is  Planck's  constant,  and  V0  is 
an  inner  potential  which  defines  the  zero  of  kinetic  energy. 

Pressing  either  a  plot  button  (for  CALCOMP  or  SC4020  plots)  or  the 
CONTINUE  button,  returns  the  menu  display.  At  this  point  select 
line  two  of  the  menu  and  type  the  upper  and  lower  transform  limits 
In  k  with  the  standard  keyboard  or  press  continue  and  CHI  versus 
k  will  again  appear  to  allow  selection  of  the  data  range  to  be  trans¬ 
formed  (using  the  light-pen)  or  choose  any  combination  of  typed  limits 
or  light-pen  are  automatically  marked  with  a  "-f".  This  display  may 
be  plotted  and  again  the  menu  returns  after  either  plotting  or  a 
continue  is  pressed.  The  initial  data  CHI  vs.  k  for  crystalline  Ge, 
with  a  typical  transform  range  selected  by  the  light-pen  is  shown  in 
Figure  2. 
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The  FFT  algorithm  requires  input  data  which  is  equi -spaced  in  k. 

This  is  generated  by  touching  the  next  line  of  the  menu  with  the 
light-pen  which  uses  ?  parabolic  interpolation  routine  to  equi-space 
the  data  (Chibar  vs.  Kbar).  The  number  of  points  displayed  in  the 
range  of  the  transform  is  proportional  to  the  number  of  input  data 
points.  These  points  are  followed  with  zeros  to  provide  a  total  of 
2^3  complex  data  pairs.  There  are  400-1000  values  of  Chibar  in  a 
typical  transform  range  and  the  remaining  7000  or  so  points  are 
zeros.  This  technique  (Bingham,  et  al.,  1967)  is  widely  used  to 
achieve  higher  resolution  in  the  transform.  By  extending  the  length 
of  the  period  in  k-space  through  adding  zeros,  the  distance  between 
harmonics  (r-space)  is  reduced.  The  transform  input  files  are  filled 
and  the  display  for  Chibar  vs.  Kbar  (the  equi-spaced  values  of  k), 
appears  on  the  screen.  The  user  has  an  option  of  multiplying  Chibar 
by  a  convergence  factor  A-exp(-B*k2)  where  default  values  are  A-l 
and  B=0.  A  plot  of  CHIBAR  vs.  KBAR  with  a  convergence  factor 
(a2*B«.01)  is  shown  in  Figure  3.  This  data  may  also  be  printed  O"  the 

line  printer. 

The  next  line  performs  the  actual  transform.  While  this  operation  is 
progressing,  the  function  buttons  and  scope  display  are  off  to  save 
core  storage.  Typical  transforms  take  about  3  minutes  depending  on 
other  multiprocessing  demand.  We  are  currently  improving  this  version 
to  reduce  this  time  to  less  than  1  minute.  When  the  transform  is 
complete  the  menu  reappears  on  the  scope  and  the  DISPLAY  TRANSFORM 
option  is  selected.  An  option  is  presented  to  multiply  the  transform 
by  CrN  (default  values  A=1,N=0).  The  real  part,  imaginary  part,  and 
magnitude  (magnitude  -  [(real)2  *  (imag)2]1/2  of  the  Fourier  transform 
between  r=0  and  10  A  are  consecutively  displayed  by  pressing  the  plot 
or  continue  buttons.  This  data  may  also  be  printed.  The  real  and 
Imaginary  part  as  well  as  the  magnitude  of  the  transform  of  the  data 
from  Figure  3  are  shown  in  Figure  4  where  Cr  *  irr2; 
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Radial  Distance,  r  (A) 

Figures  4a,  b  and  c  The  real  part,  imaginary  part 
and  magnitude  of  the  FFT  for  the  data  of  Fig.  2.  The 
vertical  scale  is  proportional  to  the  number  of  atoms 
on  each  coordination  shell  and  has  been  multiplied  by 

> P-.  The  horizontal  scale  is  the  radial  distance  from 

o 

the  absorbing  atom  in  A. 


tude  vs  radial  distance,  r. 
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After  completion  of  this  step  the  menu  is  displayed  again  on  the 
screen.  The  options  on  lines  8-10  allow  all  or  part  of  the  s- 
form  in  r-space  to  be  Fourier  transformed  back  to  k-space. 
very  useful  in  structural  analysis  and  other  applications .because 
allows  the  contribution  of  a  particular  feature  or  range  in  r-space 
studied  in  k-space  without  interference  from  other  features, 
from  line  3  of  the  menu,  the  range  of  data  to  be  inverse  transfonned 
is  selected  using  the  typewriter  or  light-pen  as  was  done  nth  he 
initial  transform.  All  data  below  the  lower  limit  in  r  a 
the  upper  limit  in  r  are  set  to  zero.  No  interpolation  is  neede 
because  the  data  from  the  original  transform  is  equi-spaced  in  r. 

Selection  of  the  next  line  performs  the  inverse  transform  yust 
was  done  previously,  and  the  next  line  of  the  menu  displays  the 
St  Figure  5  shows  the  result  of  isolating  a  peak  in  r-space  (the 
f  rst  peak  in  Figure  4CJ  and  subsequent  inverse  Fourier  transforms  ion 
This  shows  what  part  of  the  signal  Chi  vs.  k  contributed  to  the  selecte 
0  n  -spao.  A  comparison  of  Figure  5  and  Figure  2  shows  whic 
features'  in  the  original  data  can  be  identified  with  the  first  peak. 

,f  „„  other  inverse  transforms  are  required,  the  user  has 
of  taking  another  transform  over  a  different  r  9  N 

-  -  — rc;  rr;  cr.  *  -- 

ShUTd  take  less  than'lO  minutes  in  actual  time  on  the  computer 
SO  that  several  cases  per  hour  can  be  analyzed. 

He  have  found  this  program  invaluable  as  an  aid  for  Pressing  FXAFS 
data.  This  interactive  program  package  is  availa  e  or  9 
at  Boeing  Computer  Services,  Inc.,  Seattle,  M^hin^  '  in 

experience  we  have  found  certain  opera, ^  with  . 
advance  and  repetitious  in  nature.  These  problems 
similar  batch  program  which  does  not  require  user  inter.  • 
program  is  also  available. 

*~^^^^^'^°b^ethenuPsrtArmy^RosearchCoff  ic^unde/contract^DAHCOA-Tl-COOlO 

and  monitored  by  the  U.S.  Army 
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Convert 
to  Chi  Vi  K 


Display  Chi  vs.  K 
Choose  Umits  tor  Transformation 
Display  Chi  vs.  R  for  Range  to  ba  Transformed 
Perform  Transform 
Display  the  Transform 
Read  Next  Data  Casa 


User  Flow  Chart  for  Interactive  Graphics  Program  ‘ 


Not*.  Hardcopy  SC4025  or  Calcomp  Plots  Available 


for  Every  Scope  Display 


User  Selection 
Of  T ransf r>',Tl 
M'.itipi  cr  Constants 
C-r-N 

Da*  suit  A*  1  N*0 


Display 
Chi  vs.  K 
For  Data 
Just  Read 


Allow  User 
to  Type  in 
Upper  &  Lower 
Limits  in  K 
For  Transform 


^Type  iriv^ 
.Both  Limitsjs 


Display  Chi 
vs.  K  With 
Selected  or 
Default  Limits 
Ma'ked  + 


User  Selection 
of  Massage  Constants 
A-  Exp  l-3-K^l 
Oefault  A-1.J3-0 
to  Multiply  ChT 


Perform  FFT 


Store  Phase 
Corrected 
Results  on 
Aaditional  Filet 
(59  &  60) 


Display  OiT  vs.  H 


Display  Real  Vs.  R 
Multiplied  by 
C*r  Taken 
From  Files  59  &  60 


Display  Imaginary 
vt.  r  Multiplied 
By  C-r .*N  From 
Files  59  &  60 


User  Selects 
Uoper  &  Lower 
Limits  Not  Pre¬ 
viously  Typed 
in  Using  Light 
Per.  on  Display 


Display  Magnitude 
VS.  r  Multiplied 
By  C  ■  r*.N 
r rom  Files  59  &  bu 


Chi  vs.  K  ? 


Transform 


(Displays) 


Allow  User  to 
Change  Default 
Limits  of  r  -Q..  10. 
by  Type  in  From 
Keyboard 


1 

1 

2 

3 

Display  Previous  Menu 

Select  Range  for  Non-Zeroed 

Values  for  Inv.  FFT. 

• 

4 

Display  Inverse  FFT. 

2 

3 

4 

Did  User* — 
Change  .No 

Either  * 

Limit  7  / 


Display  Magnitude 
vs.  r  as  Last 
Displayed.  +’s 
Mark  Defa  ilt  Limits 
Select  One  or  Both 
Limits  With  Light 
Pan 


Perform 

Inverse 

FFT 


Writa  Zeros  on 
Unmassaged  But 
Phase  Shifted 
Transform 
Filling  Transform 
Input  Files 
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